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^ , Abstract 

Sine-Gordon kinks are a much studied integrable system that possesses 
multi-sohton solutions. Recent studies on sine-Gordon kinks with space- 
dependent square-well-type potentials have revealed interesting dynamics of 
a single kink interacting with wells and barriers. In this paper, we study a 
^ ■ class of smooth space-dependent potentials and discuss the dynamics of one 

r — I , kink in the presence of different wells. We also present values for the criti- 

^P I cal velocity for different types of barriers. Furthermore, we study two kinks 

interacting with various wells and describe interesting trajectories such as 
double-trapping, kink knock-out and double-escape. 
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1 Introduction 

Topological solitons arise in many classical field theories [18] . One of the simplest 
systems to admit solitons is the sine-Gordon model. These solitons are known as 
kinks and occur in (1+1) dimensions. The model is described by the Lagrangian 
density 

where (p{x,t) is a scalar field and A is a coupling constant. Applying the standard 
variational principle leads to the field equation 

^-^ + Asm^ = 0. (2) 

This so-called sine-Gordon equation has applications in a large number of areas of 
physical and bio-physical interest, including the Josephson effect [H], nuclear physics 
[20] . non-linear optics [I9], ferromagnetic spin chains [25] and wave propagation in 
brain microtubules [1], amongst others. Taking a spin chain in the easy plane model 
of a ferromagnet as an example, A is determined by the strength of the magnetic 
field, which occurs in a direction in the (x-y) easy plane (see [25] for details). When 
this system becomes inhomogeneous, that is the magnetic field has a 2;-dependence, 
then the corresponding parameter A is also space-dependent. The Bloch equation 
describing the system still leads to a sine-Gordon equation, and solitons can still 
occur. We will not discuss physical interpretations further. In this paper, we shall 
be interested in the sine-Gordon model with parameter A, which, as a function of 
space, is of the form of a smooth well or barrier. For our system, A will be positive 
and, far away from the origin, very close to unity. 

The interaction of a sine-Gordon soliton with a potential obstacle was first inves- 
tigated by Fei et al [5]. The authors investigated a kink incident on a point defect 
and found a novel behaviour in which, for certain incoming velocities, the kink is 
reflected backwards, due to its interaction with the defect. Many of the results in [9] 
were explained by Goodman and Haberman [12] in terms of a two-bounce resonance 
model. In [22] Piette and Zakrzewski investigated a sine-Gordon kink interacting 
with a potential with a space-dependent term in the form of a rectangular well. They 
found that back-reflection for some incoming velocities also arises for this system 
and gave two effective models to account for this phenomenon. Soliton interactions 
with rectangular wells have also been discussed for various other soliton systems 
[21 in El El [IDl [HI [21] • A kink interacting with a smooth well is discussed in [13] in 
terms of the sine-Gordon model and in [TS] for the A(/)^ model. While these papers 
address the dynamics of a kink interacting with a well, multi-kinks interacting with 
wells have so far not been investigated. 

In this paper, section [2] gives a review of sine-Gordon kinks. Then section [3] 
introduces a two-parameter family of space-dependent potentials which include a 
variety of wells and barriers. In section H] the results of simulations of the dynamics 
of a single kink interacting with a particular well are given. Following this, in section 
[5l the dynamics of a moving kink incident on a well with another initially at rest in 
the well are presented. The paper also gives plots of the scattering data for 1-kink 



and 2-kink systems with two different types of wells, these are discussed at the end 
of section O The paper ends with a conclusion. 

2 Sine-Gordon kinks 

In this section, we recall some basic facts about sine-Gordon kinks and set up our 
notation. 

The energy of static sine-Gordon kink is 

^= f (lit) +Ml-cos</>)ldx. (3) 



One can also define a super potential, W, given by 

A(l-cos(/.). (4) 



1 /dW^ ^ 



2 
The energy (|3]) can then be rewritten in terms of two integrals 



1 r fdS dW^ ^ 



-oo </'{-oo) 

For a minimum energy soliton, we obtain the Bogomolny equations 

^ = ±— (6) 

dx dcj) ' 

where the ± sign gives a kink or anti-kink. Substituting the solution of (jl]) into (E]) 
gives the kink field 

0(x) = 4 arctan f exp f VX{x — xq) j j , (7) 

where xo is a constant of integration and corresponds to the position of the kink. 
The energy can then be evaluated to be 8vA, so one has the Bogomolny bound 

E > 8v^. (8) 

Since the Lagrangian density ([T]) is Lorentz invariant, a boost can be applied to ([7]) 
to give a moving kink with field configuration 

= 4 arctan (exp (7 {x — vt — Xq))) , (9) 

where v is the velocity of the kink, with — 1 < f < 1, and 7 = ,/_ ^ is the Lorentz 
factor, and we have set A = 1 for simplicity. There is a topological charge A^ 
associated with a sine-Gordon kink 

00 

N=^ [ ^dx. (10) 

271 dx ^ ' 



At X = ±cxD the field (f) needs to minimize the potential energy in (j3]), hence we 
choose (;/)(— oo) = and (t){oo) = 2tiN . Therefore, equation flTUl) clearly takes integer 
values. For charge N the Bogomolny bound becomes 



0(oo)=27rAf 



E> I dW = %\N\^. (ir 



0(-oo)=O 

For the case of multi-kinks this bound cannot be saturated as there is a repulsive 
force between two static kinks. Indeed, in |20] Perring and Skyrme show that the 
asymptotic interaction energy of two sine-Gordon kinks is given by 

E,„t = 32v^e^^«, (12) 

where R is the separation between the kinks. 

As well as being Lorentz invariant, the sine-Gordon model is also fully integrable 
(see [23] for a review of integrable systems). Therefore, it has the feature that 
solution generating techniques, such as the Backlund transformation, are applicable 
(see e.g. [2]). For N = 2 this transformation can be used to obtain the field 

= 4 arctan , V"^ . (13) 

I cosh ^^^ ' 



This describes two individual kinks, one static at negative infinity in time and the 
other travelling with speed v at this instant. The field configurations of the kinks 
gradually approach each other with time and the solitons interact. For simplicity, 
we have set A = 1 in formula (fT3l) . 

In the sine-Gordon model there is also a breather solution, which exists in the 
charge zero sector and can be interpreted as a kink anti-kink bound state. A solution 
that is a bound state of a breather and kink, the sine-Gordon wobble, is described 
in [16]. These solutions show interesting behaviour when they are interacting with 
square wells [TDl [H] . 

3 A class of space-dependent potentials 

In the following, we describe a class of space-dependent potentials, where the cou- 
pling constant A in ([2]) becomes space-dependent. Square- well-type potentials have 
been discussed for example in [22]- These have the disadvantage that they are not 
smooth. Here, we propose a smooth two-parameter family of barriers and wells, 
with the additional advantage that there is an analytic solution for a static kink 
located at the centre of the barrier or well. Figure [1] displays the wells which we 
study in this paper, and Figure 12] shows a selection of barriers. 

The static kink ([7|) located at x = Xq solves the static field equations 

a,,0-Asin(0) = O, (14) 




Figure 1: Wells with A(0) = | for various values of a and b. Detailed scattering 
experiments have been performed on the wells displayed with solids lines. These are 



the quartic well (a = ^^, b = ^^), a narrow well without humps (a = 0, 6 = vlO/4) 
and a generic well with humps (a = —0.25,6 = 1.126166088). Some other examples 
of wells are displayed in dotted lines, for example the double well (a = |, & = |). 



for constant A. By allowing A to depend on x we can find a two parameter family 
of kinks located at x = which solves the static field equations flT^ . The kinks are 
given by 

(f){x) = TT + 2 arctan {ax + 6sinh(a;)) 

provided that we set 



where 



and 



A(x) 



n2X + riix + riQ 
d^x^ + d2x'^ + diX + do ' 



-a^6sinh(2;) 
a' 



n2 = 

rii = 2a (6^ + a^ + 2a6cosh(x)) , 

no = 6sinh(x) (6^ + 2a^ — 1 + 4a6cosh(x) + b^ cosh(a;)^) 



4 = a^ 

(^2 = 3a^6sinh(a;), 

di = a {1 - 3b'^ + 3b'^ cosh{xf) , 

do = b sinh(a;) (l - 6^ + b^ cosh(x)^) 



(15) 
(16) 



Note that the trivial vacuum solutions </)(x) = 27m also satisfy flT^ for all functions 
X{x) provided n is an integer. 

It is easy to see that X{x) = A(— x) and that the map 

(a, 6) H-> (-a, -6) (17) 

leaves A(x) invariant. Physically, this means that A is the same for a kink (6 > 0) 
and an anti-kink (6 < 0). In the following, we restrict our attention to kinks and 
set b > 0. Since the sinh(a;) term dominates for large |a;| the parameter a can take 
positive and negative values provided fe > 0. Setting b > also guarantees that 

lim A(x) = 1, (18) 

|a;|— >oo 

which means that the kink located at xq given by (I7l) tends to the exact solution for 
large \xo\, that is when the kink is far enough away from the well. For fe = 0, the 
asymptotics change. In this case, 

A(x) = ^^, (19) 

which tends to as |x| — )■ oo. Since we want the usual sine-Gordon kink ([7]) to be 
an asymptotic solution, we only consider the case b > 0. 

In order to determine the values of a and b for which A(a;) is regular, we solve 
equation ( 1T^ for A and obtain 

X{x) = ^. (20) 

Therefore, singularities can only occur for sin = 0. With our choice of 6 > we 
have 

< arctan (ax -|- osinh(x)) < — , (21) 

hence, since arctan is a monotonic function and arctan(O) = 0, A(x) can only have 
singularities at g{x) = 0, where 

g(x) = ax + bsmh.{x). (22) 

Taking the derivative of (l22l) leads to 

^'(x) = a + 6cosh(x). (23) 

This implies that g'{x) > for a + 6 > 0, so that g{x) is a monotonic function whose 
only root is x = 0. We can evaluate the limit 

, ^, , 2(a + b)^ -b 

limAx =^ '- . 24 

x^o ^ ' a + b ^ ' 

Hence, A(x) diverges for b > and a — )■ —b. For a + b < 0, the function g{x) has 
one maximum and one minimum and hence two non-trivial zeros. Generically, A(x) 
will be singular for a + b < 0. 



Physical applications usually demand that A(x) > 0. This condition is satisfied 
provided a satisfies the stronger inequality 

a + 6>Qy. (25) 

Furthermore, we are mostly interested in wells located at the origin, which implies 
A(0) < 1. 

An interesting one parameter family is obtained by setting a = 0. It follows from 

(EH) that A(0) = 262-1. Therefore, we obtain a well with \{x) > 0, for J^ < 6 < 1, 
whereas for 6 > 1 this is a barrier. This one-parameter family is an example of a 
pure well or barrier, respectively, because the only extrema of A(x) is at x = 0. 

Fixing the value of A(0) = h and h gives a cubic equation for a, which is most 
conveniently expressed in terms of the variable a + 6, namely 

(a + 6f-^(a + 6)-^ = 0, (26) 

This cubic has the discriminant 

D = + — . 27 

216 16 ^ ^ 

Hence equation (!26|) has three real solutions for D < 0, two real solutions ioT D = 
and only one real solution for D > 0. However, not all solutions will lead to regular 
wells. 

In order to gain a better understanding of the various wells and barriers it is 
useful to calculate 

,. £Xix) -12(a + 6)6 + 126(a + 6)3-6(a + 6) + 62 

lim ; — - — — ; -— . (2o) 



X 



When this limit vanishes we have a quartic well. 

In Figure m we considered wells with A(0) = |. There is a one-parameter family 
of wells with humps on both sides for a < 0. Our main example is a = —0.25 and 
b = 1.126166088. There is another one parameter family of wells with a > which 
includes pure wells, the quartic well and double wells. As can be seen in Figure [H 
the pure wells for a > are wider than the pure well at a = 0. Using equation fl28|) 
we can derive the values for the quartic well, namely, a = ^^ and b = ^^ with 

A(0) = ^. For a > ^^ we obtain double wells, where there is a local maximum at 
X = 0, and global minima at either side. 

In Figure [3l we show interesting examples of barriers. Again, there are pure 
barriers, barriers with wells at either side and volcano (or double peak) barriers. 

4 Dynamics of one kink in the presence of wells 
and barriers 

In this section, we consider a sine-Gordon kink travelling towards our standard well 
with various initial velocities and plot the trajectories. We also calculate the critical 



velocities for different barriers and compare these to an analytic approximation. 
First, we briefly comment on our numerical scheme. 

The equations of motion for the kinks have been solved using a standard fourth 
order Runge-Kutta method with gridsize 10001. Plus and minus "infinity" are 
located at 50 and —50, respectively, so that the stepsize in space is Ax = 0.01. 
The stepsize in time has usually been taken to be At = 0.0001 which is the same 
choice of parameters as in [22] • For small initial velocities of the kink, larger values 
of At are appropriate. 



X -10 




(a) V = 0.2 





150 200 250 



300 



(b) V = 0.33 




(c) V = 0.3873 



(d) V = 0.5 



Figure 2: Kink trajectories x{t) for a well with a = —0.25 and b = 1.126166088. 



Figure 2(a) shows elastic back- reflection. In Figure 2(b) the kink becomes trapped 



in the well. In Figure 2(c) the kink is back-reflected by the well. In Figure 2(d) the 
kink crosses the well. 



4.1 Plots of trajectories 

Figure [2] illustrates all relevant phenomena for one kink interacting with a well. We 
have chosen the well with parameters a = —0.25 and b = 1.126166088, which has a 
minimum A(0) = | and two humps at x ~ ±2.20 of height 1.16. 
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For low velocity v < 0.3, there is elastic back-reflection from the hump, so that 



initial and final velocity are equal in size, as in Figure 2(a) For higher velocity 



the kink overcomes the initial barrier and becomes trapped in the well which is 



illustrated in Figure 2(b) The critical velocity Uc is defined as the smallest initial 
velocity which allows the kink to cross the well (or barrier). For the well discussed 
here, the critical velocity is Uc ~ 0.387. An example of a kink crossing the well is 



given in Figure 2(d) where it can be seen that the kink gains speed inside the well. 



The most interesting behaviour happens for a narrow range of velocities just below 



Uc- As Figure 2(c) shows, the kink may enter the well, get reflected from the well 
and then leave the well travelling in the opposite direction. This type of behaviour 
has been discussed in detail in [22]. For our standard well, this resonance window 
is very narrow. In fact the kink is trapped for v = 0.38735, is back-reflected for 
V = 0.38736 and trapped once more for v = 0.38737. For v = 0.38738 the kink 
escapes. 

4.2 Critical velocities for one kink interacting with a barrier 

In this section we follow an argument in [6] to calculate an analytic approximation to 
the critical velocity, Uc- We then compare the analytic calculations with numerical 
results. 

Making use of the moduli space approximation [17], we make the ansatz 

0(x; X) = A arctan (exp (x - X(t))) , (29) 



where X{t) is the position of the solution as a function of time. Substituting (!29|) 
into ([1]) we obtain 

gj^2g2(x-x) 8e2("="^) 

^ = (i + e2(.-x))2 - (1 + ^(^)) (i + e^(-^))^ - ^^°^ 

Note that A is now a function of x. The Lagrangian is then 

oo 

L = MX-' - 1) - y (3;,,v,). d^- (31) 

— oo 

This has a non-trivial potential term, and the motion of a soliton on the moduli 
space is therefore not geodesic for non-constant A(x). This is expected as there will 
always be a static force on the solution when in the presence of a well or barrier 
such as the one in the model presented. 

The total energy of a soliton can be calculated 

oo 

E^lsr^ + i+j ^^'['^jZl7y d.. (32) 

— OO 

In the absence of an obstacle, A(a;) = 1, and the energy is 

E = -8X^ + 8. (33) 

9 



For a soliton at rest i.e. X = the kink mass is 

Mrest = 8. 



(34) 



In the case of a pure barrier, A(x) > 1, therefore, for a kink stationary on top of the 
barrier 

M,e.t > 8. (35) 

For a pure well, on the other hand, A (a;) < 1, so a static kink at the bottom of the 
well will have mass 

Mrest < 8. (36) 

Otherwise, one cannot say which inequality holds without explicitly evaluating the 
integral appearing in (152]) . Energy conservation dictates that in the case of elastic 
scattering the mass, Mb, of a static kink on top of a barrier, is related to the critical 
velocity, Uc, for the barrier by the following formula (appearing in [5]) 



VT 



M 



B- 



Ur 



(37) 



Figure [3] shows some barriers of different types and heights. 




Figure 3: Different types of barriers: a barrier with two side wells (a = 0.25, h = 
0.75), a volcano barrier (a = —0.5, b = 1.657970214), and three pure barriers (a = 0, 
b = ^), (a = 0, 6 = ^) and (a = 0, 6 = ^2) with height 1.25, 2 and 3, respectively. 
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In table [T] we present critical velocities calculated both numerically and using 
equation fl57|) . It is clear that the theoretical approach is accurate for barriers 
that are monotonically increasing before the origin and monotonically decreasing 
afterwards. Otherwise, there is discrepancy. This is explained by considering the 
fact that for all models investigated so far solitons have elastic scattering on barriers 
and inelastic scattering on holes. 



a 


h 


mass 


Uc (theoretical) 


Uc (numerical) 


A(0) 


Type of barrier 


0.25 


0.75 


8.33 


0.271 


0.280 


1.25 


Barrier with wells 


-0.50 


1.66 


9.17 


0.489 


0.505 


1.25 


Volcano 


0.00 


1.06 


8.64 


0.379 


0.379 


1.25 


Pure barrier 


0.00 


1.22 


10.35 


0.635 


0.636 


2.00 


Pure barrier 


0.00 


1.41 


12.28 


0.759 


0.762 


3.00 


Pure barrier 



Table 1: The values of a and h are the parameters of the barriers discussed in Figure 
[31 the column "mass" gives the energy of the a kink located on top of the barrier, 
the two columns with Uc compare the analytic approximation to the critical velocity 
to its numerical value, A(0) is the height of the barrier at x = and the final column 
gives the type of barrier. 



5 Dynamics of two kinks in the presence of a well 

In general, the interaction of two kinks with a well is rather complicated because 
the trajectories depend on the initial positions and the initial velocities of both 
kinks. Here, we restrict our attention to scattering processes when the first kink 
approaches the well with given initial velocity w, and the second kink is at rest 
in the well. Our potential A(x) has been chosen, such that the first kink is given 
asymptotically by equation (|9]) whereas the second kink is given by equation f lTSj) . 
The exact solution flT^ cannot be generalized to our wells. However, since kinks 
are exponentially localized, we can just concatenate the two kinks provided the first 
kink is far enough away from the well. 

Figure H] illustrates some important trajectories for our standard well. When 
the first kink travels with small initial velocity, e.g. v = 0.2, then there is elastic 
back-reflection, such that the final velocity of the first kink is —v, and the second 
kink remains at rest in the well, see Figure 4(a) The trajectory of the first kink 



closely resembles the trajectory of a single kink in Figure 2(a) Figure 4(b) shows 



inelastic scattering where the first kink loses energy, and the second kink clearly 
gains energy because is it now oscillating in the well. In Figure 4(c) the second 



kink is knocked out of the well by the first kink which also escapes from the well. 



Figure 4(d) shows a similar scenario. The second kink is again knocked out of the 
well, whereas the first kink is now trapped in the well. It is worth comparing this 
outcome to the exact solution when there is no well. In this case, the first kink loses 
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(b) V = 0.6 




(c) V = 0.7 



(d) V = 0.75 



Figure 4: Kink trajectories x{t) for a well with a = —0.25 and b = 1.126166088. One 
kink is travelling towards the well whereas a second kink is stationary inside the well. 
Figure 4(a) shows elastic back-reflection, the kink in the well remains stationary. In 



Figure 4(b) the kink experiences inelastic scattering, the trapped kink oscillates in 
the well. In Figure 4(c) the flrst kink is back-reflected while the second kink escapes 



from the well. In Figure 4(d) the flrst kink becomes trapped while the second kink 
escapes from the well. 



all its kinetic energy to the second kink, so that after the scattering event the flrst 
kink is at rest whereas the second kink moves with velocity v. 

We checked the range 0.695 < v < 0.6951 in detail, in the hope that there would 
be a double back-reflection, so that flrst kink and second kink both escape on the 
same side. But for our standard well, we were unable to flnd this kind of behaviour. 

In Figure [5l we discuss two interesting trajectories for the quartic well (a = ■^^, 



32 



This is a rather wide well which allows both kinks to be trapped. Figure 
5(a) shows an interesting quasi-periodic motion where the two kinks knock each 
other out of the centre of the well but not gaining enough energy to leave the well. 



Figure 5(b) shows a novel type of back-reflection. The flrst kink enters the well and 
speeds up, then hits the second kink and knocks it off the centre of the well. The 
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Figure 5: Kink trajectories x{t) for the quartic well (a = ^^, b = ^^)- One kink 
is travelling towards the well with velocity v whereas a second kink is stationary 



inside the well. In Figure 5(a) both kinks become trapped in the well. Figure 5(b) 
shows a novel type of back-reflection. 



second kink slows down, returns to the centre and gives the first kink enough kinetic 
energy to escape the well. 
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Figure 6: Initial velocity versus final velocity for a kink incident on a well with 
humps with a = —0.25 and b = 1.126166088. For the figure on the left, the well 
is empty whereas for the figure on the right, there is a kink in the well, which is 
initially at rest. 



Figure [6] shows initial versus final velocity for our standard well. In the left 
figure, one kink travels towards the well with initial velocity v. For v < 0.3 the kink 
is elastically back-reflected. As v is increased, the kink is trapped. For v > 0.387 the 
kink escapes the well. As v tends to one the kink feels the influence of the well less, 
so that the final velocity tends to the initial velocity for v ^ 1. The back- reflection 
at V = 0.38736 is too narrow to be detected with the resolution used for this figure. 
However, there is another backreflection at v = 0.3860. 
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The right figure in|6]shows initial versus final velocity for our standard well, where 
a second kink is at rest inside the well. Elastic scattering takes place for a longer 
range of initial velocity v. As the initial velocity of the incoming kink increases, the 
kink in the well is set oscillating more energetically. Because of this, at some point 
the final velocity of the first kink starts to become less negative. At higher initial 
velocities the final velocity of the initially static kink increases and approaches one, 
while the final velocity of the other soliton takes the value zero in two separated 
"windows" , but is more negative between and after these windows. 

Figure [7] shows initial velocity versus final velocity for a narrow well without 
humps. In the left figure, one kink travels towards the well with initial velocity 
V. There is a "resonant window" at v = 0.22 where back-refiection occurs. In 
comparison to Figure E] the behaviour near the critical velocity is less abrupt. For 
large initial speed, the initial velocity approaches the final velocity. However, there is 
a region where the curve appears to "wobble" . In terms of our numerical simulations, 
we have checked sensitivity to stepsize in time and space independently. The wobble 
is more sensitive to the stepsize in space, but appears to be a genuine phenomenon. 
This could be a novel feature due to our potential. 

The figure on the right in [7] shows initial velocity versus final velocity for the 
narrow well without humps, where a second kink is initially static in the well. At 
small to medium initial velocities, elastic and in-elastic back-refiection occurs. For 
velocity v > 0.6 the static kink is ejected from the well and at marginally greater 
initial velocity, the first kink becomes trapped. For higher velocities both kinks 
escape the well in opposite directions. Again the final velocity of the first kink is 
not a monotonic function of the initial velocity. 
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Figure 7: Initial velocity versus final velocity for a kink incident on a narrow well 
without humps with a = and b = vTO/4. For the figure on the left, the well 
is empty whereas for the figure on the right, there is a kink in the well, which is 
initially at rest. 
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6 Conclusion 

In this paper, we have discussed the interactions of kinks with wells and barriers. 
We introduced a novel class of smooth potentials which included pure wells and 
pure barriers of various heights and widths. These potentials also give rise to wells 
with two side barriers and barriers with two side wells, as well as double wells and 
volcano barriers. The main advantage of the proposed potentials is that the static 
kink-in-the-well solution is explicitly known. The asymptotic kink solution for a 
kink located at very large xq is also known. Therefore, a kink can be scattered off a 
kink-in-the-well and the only relevant parameter is initial velocity of the incoming 
kink. When the exact solution of the kink in the well is not known, then the trapped 
kink tends to oscillate in the well, such that the scattering behaviour depends not 
only on the velocity of the incoming kink but also on the phase of the trapped kink, 
which is more difficult to control. 

We studied the scattering of one kink off a well, mainly focussing on our standard 
well with A(0) = | and two side humps, and reproduced all the known phenomena 
such as trapping, back-reflection and escape. We then compared numerically calcu- 
lated critical velocities for various barriers to an analytic approximation assuming 
elastic scattering. As expected there was good agreement for pure barriers, and less 
agreement for the volcano barrier and the barrier with side wells. It would also be 
interesting to consider the two-kink system, with one kink initially static in a well 
and to compute the "critical velocity" of an incoming kink at which this static kink 
is knocked out of the well. 

The novel feature of this paper is the scattering of multi-kinks in the presence 
of a well. We calculated various trajectories for our standard well, such as elastic 
scattering, inelastic scattering, kink one replacing kink two in the well and a scat- 
tering where both kinks escape. Unfortunately, we were unable to find a double 
back-refiection for our standard well, but it is likely that such a trajectory exists 
for the right choice of well parameters and initial velocity. For the quartic well, we 
found double trapping of kinks and a novel type of back-refiection. We also plotted 
initial velocity against final velocity, both for an empty well and a kink trapped in a 
well. In the case of our standard well with humps, we found that, for a single kink, 
the plot reproduced all the expected behaviour. In the case of an incoming kink and 
a kink initially in the well we found for small to medium initial velocities there is 
back-refiection that starts of as elastic then becomes inelastic. When the initial ve- 
locity becomes high, the static kink is ejected and there are two "windows" in which 
the initially moving kink is trapped. For a narrow well without humps the plot for a 
single kink revealed a "resonant window" and an interesting "wobbling" behaviour. 
The plots for two kinks showed that as the initial velocity increases, there is a point 
where the static kink is ejected and further on, the first kink becomes trapped. At 
high speeds both kinks escape the well, moving in opposite directions. Again, the 
final velocity of the first kink is not a monotonic function of the initial velocity. 

A long term aim of this line of work is to be able to design wells which have 
desirable properties. In certain physical systems, it is possible to make the coupling 
constant inhomogeneous, e.g. the ferromagnetic spin chain described in the intro- 
duction [25]. Barriers and wells might then be used to control kink dynamics for 
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example by acting as filters. An even more ambitious aim is discussed in [8] where 
a particular kind of integrable point defect could be used to construct simple logical 
gates. 

Our approach is applicable to various different systems. The A0^ kink allows 
similar wells and barriers which would enable the study of kink anti-kink scattering in 
the presence of barriers and wells. Soliton-well interaction has also been discussed for 
various other solitons for example deformed sine-Gordon models [U [7], A0^ models 
[3], Q-ball systems [3] or generalized sigma models [TU] where our approach may 
again prove to be useful. 

Reference |23] presents calculations for two interacting kinks using the method 
of Manton for constant A. In |22] a moduli space approximation was proposed which 
also takes account of the degrees of freedom of the well. See [T3] for related work on 
the X(j)^ kink. Our calculations imply that the an additional degree of freedom could 
be the slope of the soliton at the centre of the kink as parametrized by a and b. One 
could therefore investigate moduli space dynamics applied to the model described 
in this paper, where a, b and xq for each kink become functions of t, and the degree 
of freedom of the well is also taken into account. This is currently work in progress. 
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